IUHET 278 
May, 1994 



Predicting the masses of heavy hadrons without an explicit Hamiltonian 



R. Roncaglia, A. Dzierba, D. B. Lichtenberg, and E. Predazzi,* 
Department of Physics, Indiana University, Bloomington, Indiana, 47405 



There are striking regularities in the masses and mass differences of known hadrons. Some 
of these regularities can be understood from known general properties of the interactions of 
quarks without a need to specify the explicit form of the Hamiltonian. The Fcynman-Hcllmann 
theorem is one of the tools providing this understanding. If the mass regularities are exploited, 
predictions can be made of the masses of as yet undiscovered hadrons. In particular, it is found 
that the mass of the B* is 6320 ± 20 MeV. Predictions concerning i) excited vector mesons, 
ii) pseudoscalar mesons, iii) P-wave mesons, and iv) ground-state spin 1/2 and 3/2 baryons 
are also made. 

PACS numbers: 12.15.Ff, 12.40.Yx, 14.20.-c, 14.40.-n 



* On leave from the University of Torino, Italy 



2 



I. INTRODUCTION 

Thus far, the quark potential model has been the most successful tool enabling physicists to 
calculate the masses of normal mesons and baryons containing heavy quarks. We call attention 
to several reviews on the subject [1-8]. However, potential models suffer from the fact that, 
although motivated from QCD, so far, they cannot be derived from that theory. In this paper we 
make predictions about hadron masses with complementary methods which use general properties 
of the potential (or, more generally, of the interaction) but not its specific form. Among the 
complementary methods are: 

(1) Obtaining constraints on hadron (and quark) masses [1,9,10] from the Feynman-Hellmann 
theorem [11,12]; 

(2) Using theorems which relate the ordering of bound-state energy levels to certain properties of 
potentials [13]; and 

(3) Taking advantage of regularities in known hadron masses to obtain estimates of as yet undis- 
covered hadrons using either interpolation [9] or semi-empirical mass formulas [14,15]. 

We exploit these methods to obtain constraints on quark and hadron masses. We also provide 
new theoretical justification for the methods we use and make predictions for the masses of as yet 
undiscovered mesons and baryons. 

We devote considerable effort to making what we believe is a good prediction for the mass 
of the B* vector meson, using the Feynman-Hellmann theorem. We also discuss in some detail 
the ground-state pseudoscalar mesons and the ground-state baryons (both spin 3/2 and spin 1/2), 
because additional issues arise in these cases. We discuss only briefly the excited vector mesons 
and P-wave mesons (tensors, axial vectors, and scalars), although we give some predictions in these 
cases as well. 

II. THE FEYNMAN-HELLMANN THEOREM 

Some years ago, Feynman [11] and Hellmann [12] independently showed that if a Hamiltonian 
H depends on a parameter A, then the bound-state energy eigenvalues E(X) vary with A according 
to the formula 

8E/8X= (dH/dX), (1) 

where the expectation value is taken with respect to the normalized eigenfunction belonging to 
E. The Feynman-Hellmann theorem was applied to quarkonium physics [1,10], with A = fi, the 
reduced mass of the system. For example, Quigg and Rosner [1] applied the theorem to the 
nonrelativistic Hamiltonian H = p 2 /(2/.i) + V, where V is an interaction which is assumed to be 
flavor independent, and therefore independent of fi. Then 

dE/dv=-(p 2 )/(2n 2 ) <0, (2) 

i.e., E decreases monotonically as n increases because p 2 is a positive definite operator. Of course, 
if V depends on /i but (dV/dfi) < 0, then dE/dp < still remains valid. 

Even in the case of some many-body Hamiltonians with relativistic kinematics, the Feynman- 
Hellmann theorem may be applied to give useful information about how eigenenergies change when 
constituent masses m, change [16]. As an example, we consider a Hamiltonian H, given by 

H = EK p * + m *) 1/2 ~ + F ( ri ' ••• r ™ ; mi > ■•■ m ")' ( 3 ) 

i 

where we have let the interaction V depend explicitly on the m, (i = 1, 2...n). Taking the partial 
derivative with respect to m; and using (1), we obtain 

dE/dm, = K/(p 4 2 + m 2 ) 1 ' 2 ) - 1 + (dV/dm,). (4) 



We can see from Eq. (4) that if 

(dV/drrii) < 0, 



(5) 
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then, 

dE/dmi<0, TOj = l,2...n. (6) 

These inequalities are a generalization of (2). 

In the remainder of this section, we restrict ourselves to the case in which (5) and (6) both 
hold, so that an increase in one or more rrn leads to a decrease in E. We generalize the definition 
of to be 

M _1 = E m . rl - ( 7 ) 

We now note that an increase in one or more results in an increase in /j, as well as a decrease 
in E. Under these circumstances, a change 5/j, results in a change 5E in the opposite direction. 
Therefore, E will be monotonically decreasing with increasing /i, provided the increase in /x arises 
from an increase in one or more m*. We can state this result in the form 

SE/Sfi < 0, (8) 

if the changes in the mi are all in the same direction. The inequality (8) turns out to be a powerful 
tool for obtaining constraints on quark and hadron masses. As we shall see in Sec. V, this inequality 
holds empirically for vector mesons even in the absence of the restriction that all mj change in the 
same direction. 

For a two-body Hamiltonian of the form (3) with a flavor-independent potential, we can say 
something more. Consider £ as a function of /z and M., where M. = mi + m 2 . We can then show 
explicitly that 

dE/d/j, < 0, dE/dM < 0. (9) 

It follows that even if mi and mi change in opposite directions, if ji increases and M either 
remains constant or increases, then E decreases. This result holds even in the presence of some 
flavor-dependent interactions, including a colormagnetic interaction for vector mesons, as we shall 
see in the next section. 

III. APPLICATION TO MESON AND BARYON EIGENENERGIES 

We adopt a constituent quark picture, assuming that a meson is composed of a quark and an 
antiquark, and a baryon is composed of three quarks. We confine ourselves to hadrons containing 
u, d, s, c, and b quarks. Furthermore, we neglect any violation of isospin and let m u = m d = m q . 
As usual, we assume that the quark masses satisfy the inequalities 

m q < m s < m c < nib- (10) 

In order to apply (6) to hadrons, we need to discuss for which hadrons (5) is likely to hold. The 
interaction V can be written as Vq + V, where Vq is independent of quark flavors and V' depends 
on flavor. The term Vq is the static quark-antiquark potential, which is commonly assumed [17] 
to contain a Coulomb-like term, an approximately linear confining term, and a constant term, 
all independent of flavor. In the Fermi-Breit approximation, V contains both spin-dependent 
and spin-independent terms which are explicitly functions of flavor through the quark masses. 
However, most phenomenological treatments of quarkonia have not needed the Fermi-Breit spin- 
independent term [5], and we neglect it here. In states with zero orbital angular momentum, the 
expectation values of the tensor and spin-orbit terms of the Fermi-Breit interaction vanish, leaving 
the colormagnetic interaction as the only spin-dependent term. We write the colormagnetic term 
V cm in the form [18] 

Vcm = ~E Ai ' / W( r u) cr i -o-j/imimj), (11) 

i<j 

where &i, Oj are Pauli spin matrices, Xi, Xj are Gell-Mann SU(3) matrices, and /(r^-) (r^ = [r— |) 
are positive definite operators. 
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If V = Vq + V cm , we obtain 

(dV/drm) = E</( r «))( A » • A ^ ■ ^)/K 2 ^)- ( 12 ) 

The quantity (Aj • Aj) is negative for a quark-antiquark pair in a meson (—16/3) and for all quark 
pairs in a baryon (—8/3). Also, (<r, • Uj) is —3 if two quarks (or a quark and antiquark) are in a 
spin-zero state, and 1 if they are in a spin-one state. Now in vector mesons and spin-3/2 baryons, 
we have (ct, • aj) = 1 for all quark pairs. Then we see from (12) that the eigenenergies of these 
hadrons satisfy (5). We therefore expect the energy eigenvalues of vector mesons and spin-3/2 
baryons to satisfy (5), and therefore (6) as well. In the two-body case (vector mesons), (11) can 
be written 

V cm = (16/3)/(n 2 )/(/iA<). (13) 

We can see explicitly from (13) that (dV/dfi) < 0, (dV/dM) < 0, so that (9) holds, as we stated 
in the previous section. 

For pseudoscalar mesons the sign of the colormagnetic term is negative, and the interaction 
violates (5). For spin 1/2 baryons, the three colormagnetic terms in (11) are either all negative or 
one term is positive and two are negative, so that (5) is sometimes violated. Thus, we expect that 
pseudoscalar mesons and spin- 1/2 baryons might violate (6) for small m„ where the contribution 
from (12) is large and positive. 

In addition to V cm , terms arising from instantons may contribute to V. These terms are 
apparently important in states in which two quarks (or a quark and an antiquark) have spin and 
orbital angular momentum zero (pseudoscalar mesons and spin-1/2 baryons) [19]. Instantons tend 
to mix the wave functions of certain mesons, like the n and r\' (which contain both qq and ss in 
their wave functions, and perhaps some glueball admixture as well). Such states are unsuitable 
for our scheme, as, in order to compute the reduced mass of a system, we must know its quark 
content. 

If we confine ourselves to mesons containing q(= u,d), s, c, and b quarks, then 10 different 
ground-state vector mesons and 20 different ground-state baryons of spin 3/2 exist. Of the former, 
9 are experimentally known; of the latter, only 4 are known, none of which contains a heavy quark. 

Let a particular energy eigenvalue of a meson containing a quark i and an antiquark j be E^ , 
where we suppress an index labeling which eigenvalue we are referring to. Likewise, we denote 
an energy eigenvalue of a baryon by E%jk. If we replace a single quark by a heavier quark, then, 
because of (6), the eigenenergy of the new hadron will be smaller than that of the old one. By 
continuing this process, we obtain chains of inequalities among the Eij and also among the E^. 
We consider the longest chains of inequalities for which (6) holds. A longest chain arises when we 
start with the lightest hadron (containing only g-type quarks) and replace each of its quarks one 
at a time by the next heavier quark, each time obtaining a different hadron (i.e., we replace q by 
s, s by c, and c by b). A longest meson chain contains 7 eigenenergies, and a longest baryon chain 
contains 10. 

There are five different longest chains for mesons, all of which have the same first two and last 
two eigenenergies. At least one of the three intermediate eigenenergies differs in each of the five 
chains. One longest chain for mesons is 

Ebb < Eb c < E cc < E cs < E ss < E sq < E qq . (14) 

The four other longest chains have intermediate eigenenergies 

E C c^E CS: E C q 1 Ebs : E CS: E SSl Eb s ^E CS: E C q 1 Eb S: Ebq : E C q. (15) 

It follows that (6) holds for the mesons with eigenenergies as in (14) or (15), whereas this may not 
necessarily be the case when one quark mass increases and the other decreases. 
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We next turn to the spin-3/2 baryons, and denote their energy eigenvalues by E^. From our 
discussion we can write down a number of longest chains. One such chain is 

Ebbb < Etbc < Ebcc < Eb C s < Ebss < E css < E sss < E ssq < E sqq < E qqq . (16). 

All longest distinct baryon chains, of which we have counted 42, contain the same first two and 
last two eigenenergies. 

IV. CONSTRAINTS ON QUARK MASS DIFFERENCES 

The energy eigenvalues E^ and E^k do not include the quark rest energies. Therefore, the 
mass M of a hadron is given by 

M = E + ^TO i7 (17) 

where we have suppressed the indices on E and M. As has been pointed out by a number of 
authors [1,2,20,21] Eq. (17) can be used together with (8) to obtain constraints on quark mass 
differences in the form of inequalities. What is new in our treatment is our justification of the 
use of (8) for longest chains of vector meson and spin-3/2 baryon eigenenergies. In particular, 
this means that we do not need to use as input spin-averaged hadron masses, which are normally 
calculated in a perturbative approximation. 

In order to obtain inequalities among quark masses, we use the experimental values of hadron 
masses, including the p, K* , <f>, D*, D*, B*, B* sl J/ip, and T mesons (of the 10 vector mesons, 
only the B* is missing) and the A, £*, 5*, and ft baryons (only 4 of the 20 spin-3/2 baryons are 
known). When the isospin of a state is greater than zero, we average over the members of the 
isospin multiplet. For the qq state, we believe it is better to choose the p than the w because the 
latter might have a small admixture of ss, whereas such an admixture in the p violates isospin. 
Note that in our scheme the uo is degenerate with the p, although it is actually 15 MeV heavier. 
In fact, in not distinguishing between u and d quarks, we consider the uj and p together as just a 
single one of the 10 distinct vector mesons in our scheme, which is concerned only with masses. 
We have given an a priori reason for choosing the p instead of the u>. In the next paragraph we 
give an a posteriori reason. However, if instead of choosing the p, we choose the u, our predictions 
are not appreciably affected. All the meson masses are taken from the Particle Data Group [22], 
except the mass of the B*, which comes from two recent measurements [23,24] of the mass of the 
B s plus a measurement of m(B*) — m(B s ) (which needs confirmation) quoted in [22]. 

Using the masses of the observed ground-state vector mesons, we obtain the following inequal- 
ities: 

m s - m q > M(K*) - M(p) = 126 ± 4 MeV, 

m c -m s > M(D*) - M(K*) = 1115 ± 4 MeV, (18) 

m b -m c > M(B*) - M(D*) = 3316 ± 6 MeV. 

The errors are partly experimental and partly due to our assumption of isospin invariance. If we 
substitute the u for the p, the right hand side of the first of these inequalities becomes smaller by 
15 MeV. However, we can regain the stronger inequality (larger right-hand side) by considering 

m s - m q > M{4>) - M(K*) = 125 ± 4 MeV. (19) 

This fact gives an a posteriori justification of our decision to use the p, rather than the u>, in our 
set of vector mesons. 

One can obtain apparently stronger inequalities by using the spin-averaged values for vector 
and pseudoscalar mesons. The results in MeV are [20] 



to s - m g > 184 ± 4, m c - m s > 1180 ± 4, m b - m c > 3343 ± 4. 



(20) 
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However, the spin averaging process relies on a perturbative treatment of the spin-dependent 
forces, which may not be justified for light mesons. In fact, we shall see in Sec. VI that when 
we vary the quark masses to obtain a best fit to the data, the quark mass differences violate (20) 
but satisfy (18). Still other authors [21] have obtained inequalities among quark masses using 
a variety of assumptions, but not if the Hamiltonian contains both relativistic kinematics and a 
flavor-dependent interaction. 

As we have remarked, only 4 ground-state baryons of spin 3/2 have been observed thus far. 
These lead only to inequalities for m s — m q . The strongest of these is 



This inequality is stronger than the corresponding inequality we obtained from mesons. As we 
shall see in Sec. VIII, in order to obtain a best fit to the baryon data, we must use quark mass 
differences for baryons which are up to 35 MeV larger than the corresponding quark mass differences 
for mesons. 



Many different sets of quark constituent masses are used in the literature, most of them ob- 
tained from fits to spectroscopic data. We show in Table I a selection of these sets [7,9,14,15,17,18, 
21,25-28], including in the first row the set that we use in this work for vector mesons. (In the 
next section, we explain how we arrive at this mass set.) We can see from Table I that the sets 
of masses given in the first seven rows satisfy the inequalities given in (18), whereas the sets of 
masses given in the last five rows violate one or more of these inequalities. Therefore, if one cal- 
culates the vector meson masses with any one of the last five sets and a potential which is flavor 
independent except for a conventional colormagnetic term, one will obtain results in disagreement 
with experiment. This disagreement will occur independently of whether one uses a nonrelativistic 
Schrodinger equation or a wave equation with relativistic kinematics of the form given in Eq. (3). 

We use the experimental values of the vector meson masses together with the input values 
of the quark masses to calculate the meson eigenenergies Eij from Eq. (17). We can plot these 
eigenenergies as a function of /x. For illustrative purposes, we use four different sets of input quark 
masses from Table I, and show in Fig. 1 how E varies as a function of \i for the ground-state vector 
mesons. For completeness, we include mesons containing qb, qc, and sb as well as sc., ss, and cc 
in order to see how E varies with \i when mi increases and mi decreases. In Fig. 1 (a) and (b), 
the quark masses satisfy all the inequalities (18), whereas in Fig. 1 (c) and (d), the quark masses 
violate at least one of these inequalities. 

We see from Fig. 1 (a) and (b) that, when the quark masses satisfy (18), E appears indeed 
to be a monotonically decreasing function of /i for the vector mesons. From Fig. 1 (c) and (d) we 
see that, with quark masses violating (18), E, as obtained from the observed meson masses, is not 
monotonically decreasing as a function of zz. The sets of quark masses used in (c) and (d) seem a 
priori as reasonable as those of sets (a) and (b). However, if we use the quark masses of (c) or (d) 
with a flavor-independent potential and a conventional colormagnetic term, the calculated values 
of E must be monotonically decreasing and therefore in disagreement with experiment. 

We see from Fig. 1 (a) and (b) that, not only is E monotonically decreasing as a function of 
ji, but is also concave upward. We cannot prove concavity from the Feynman Hellmann theorem. 
However, in the two-body nonrelativistic case, we can show for a power-law potential of the form 



m s -m q > M(S*) - M (A) = 153 ± 4 MeV. 



(21) 



V. SATISFYING THE QUARK MASS CONSTRAINTS 



V = ar , a(3 > 0, > -2 



(22) 



that the concavity condition 



d 2 E/d/i 2 > 



(23) 



is true. This result follows directly from the scaling property of the Schrodinger equation with a 
power-law potential [1]. The fact that (23) holds for a power-law potential is relevant because, as 
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has been emphasized by Martin [29,30], the quarkonium static potential can be well approximated 
by (22). If we include a Fermi-Breit term in the potential of form given by Eq. (12), (23) remains 
valid. Of course, we are interested in the curvature of E as a function of fi without concerning 
ourselves with how M. varies. In the presence of the colormagnetic interaction (11), we can say 
that the curvature is positive, provided that as \i increases, M does not decrease. It has been 
shown [31, 32] that in the two-body nonrelativistic case with a flavor independent potential, 

d 2 E/d\ 2 < 0, (24) 

where A = l//i. Not only is E concave upward as a function of /j, in a class of nonrelativistic 
two-body models, but it turns out that 

S 2 E/5n 2 > (25) 

holds empirically for vector mesons (see Fig. 1 a and b) and also for spin-3/2 baryons. 

VI. PREDICTING THE B* MASS 

Of the 10 vector mesons in our scheme, only the B* has not yet been seen. We can use 
the inequalities (8) and (25) together with the experimental vector meson masses and Eq. (17) 
to estimate the mass of the B* by interpolation. We do this by assuming that E{y) can be 
approximated by a simple curve containing only a few parameters. The simplest curve (containing 
only 2 parameters) is a straight line, but a straight line violates the concavity condition (25). We 
therefore approximate E(/j,) with a three-parameter curve. We emphasize that the functional form 
of the curve has no theoretical significance other than that it satisfies the inequalities (8) and (25). 
In addition to the three parameters of the curve, we have four additional parameters, namely, the 
quark masses m qi m s , m c , and m^. Using a given three-parameter curve, we vary the 7 parameters 
in order to obtain a best fit to the eigenenergies Eij . We then obtain the value of E\, c from the 
fitted curve and use the fitted masses m c and m& to obtain M(B*). 

We at once encounter a difficulty in our scheme: namely, that a longest chain of meson 
eigenenergies contains only 7 members, and one of these (Eb c ) is unknown. Therefore, if we use a 
longest chain, we have 7 parameters and only 6 data points so that the parameters are not uniquely 
determined. We overcome this difficulty by using all 9 known meson eigenenergies, assuming that 
(8) and (25) hold even in this case. We find that our assumption is consistent with experiment: 
namely, we can find a set of quark masses such that the meson eigenenergies, as calculated from the 
experimental values of the meson masses with the aid of (17), satisfy (8) and (25). The procedure 
of including all nine known meson masses constrains the parameters much more than using only 
six masses. However, even when we use 9 data points, it turns out that the quark mass differences 
are much more constrained than the quark masses themselves. 

Kwong and Rosner [9] previously used the interpolation method with two three-parameter 
curves (quadratic and Pade), although without any theoretical justification of the inequalities (9) 
and (25). We have used three different three-parameter curves, an exponential 

E = acxp(-n/b) - c, (26) 

a quadratic 

E = a + bfi + cfi 2 , (27) 

and a hyperbolic (or Pade) 

E = a/{n + b)-c, (28) 

where o, b, and c are parameters to be varied. In principle, these parameters are functions of the 
sum of the quark masses as well as a function of fi, but from our previous discussion, we expect E 
to be decreasing when plotted as a function of \x with a, b and c independent of any other masses in 
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the problem. We obtain comparable fits to the data with all three curves, and the meson energies 
are quite stable to our choice of functions. 

We show in Fig. 2 our fit to the vector meson energies, with an exponential curve as an 
example, and the set of quark masses given in the first row of Table I. These quark masses are 
rounded to the nearest 10 MeV and are based on a somewhat arbitrary choice of 300 MeV for the 
mass of the u and d quarks. We can get comparable fits to the data for quark masses which differ 
from our choices by 100 MeV or more, but the mass differences are much more constrained. With 
the exponential fit, the values of the parameters with our quark masses are 

a = 754 MeV, b = 1375 MeV, c = 506 MeV. (29) 

Our result for the B* mass is 

M(B*) = 6320 ± 20 MeV. (30) 

We have estimated the theoretical error partly from the spread in values obtained using the different 
functional forms for E(p) given in Eqs. (26-28) and partly from values obtained with different 
longest chains and various constraints on the quark masses. Our quoted errors in Eq. (20) and our 
subsequent predictions reflect the stability inherent in our procedures. Our predicted value of the 
B* mass is given in Table II in the first row of column 4. 

Our value of M(B*) is more stable to the choice of curve than the result of Kwong and Rosner 
[9]: 6284 < M(B*) < 6349 MeV. Perhaps one reason for this is that we differ from those authors 
in the choice of the function % 2 to be minimized. We choose 

X 2 = ]>>( M )-£(cxp)] 2 /(AM) 2 , (31) 

where E(fj,) is obtained from one of the three curves, £7(exp) are the experimental eigenenergies 
obtained with the help of (17), and AM are the experimental errors in the meson masses, except 
that we have taken a minimum error of 1 MeV and increased some errors to take isospin mass 
splittings into account. 

Other authors, using potential models, have obtained similar values of the B* mass. For 
example, Martin [30], using a power law potential, obtained a value of 6318 MeV, and Eichten 
and Quigg [33], using various potentials, obtained values between 6319 and 6343 MeV. Bagan et 
al. [34] have averaged a variety of other people's results to obtain 6330 ± 20. On the other hand, 
some authors have obtained quite different values of the mass of the B*. For example, Jain and 
Munczek [35] find M(B*) = 6277 MeV. 

VII. OTHER MESON MASS PREDICTIONS 

As we have already remarked, we do not expect the pseudoscalar meson energy eigenvalues to 
be monotonically decreasing as a function of /x. However, once we have a prediction for the B* mass, 
we can obtain estimates for pseudoscalar meson masses in other ways. We use semi-empirical mass 
formulas [14,15] for the splitting between vector and pseudoscalar states. These semi-empirical 
formulas are based on the colormagnetic interaction as given by the Fermi-Breit theory. However, 
the formulas take into account the fact that the strong-interaction coupling constant runs. The 
formulas also make empirical corrections which depend on quark masses so as to get improved 
agreement with known data compared to the Fermi-Breit formula. The semi-empirical formulas 
then may be used to predict mass splittings in cases where experimental data are absent. 

Following [15], we take the hyperfme splitting in mesons to be given by: 

M v -M s =pa s (2fj,)fj,y(m 1 +m 2 ), (32) 

where mi and m 2 are the constituent quark masses, and we have determined the two free param- 
eters p and q from a fit to the experimental splittings to be q — 0.642 and p — 1.917 (GeV) 2 ~ 9 . 
The running coupling constant is given by 

a 'W = to+vmj*rt' (33) 
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where 

t = ln(Q 2 /A 2 QCD ), 0o = 11 - 2n//3, fa = 102 - 38n//3, (34) 

with Aqcd = 100 MeV, and ray =4. We use the quark masses of the first row of Table I rather 
than the masses in [15], so that our values of the parameters q and p differ a little from those in 
[15]. 

We show in Fig. 3 the eigenenergies of the pseudoscalar mesons, where we use the data [22] on 
pseudoscalar masses as input together with the quark masses obtained for the vector mesons. We 
omit the rj and rj mesons, because, as we have already remarked, they are mixed states of uncertain 
quark content. We also give in Fig. 3 the eigenenergies of the pseudoscalars obtained from the 
corresponding vectors with the aid of the semi-empirical mass formula (32). We see from Fig. 3 
that, except (as expected) for the pion, the eigenenergies obtained from the observed masses and 
from the semi-empirical formulas are in remarkably good agreement. Note that the eigenenergies of 
the light pseudoscalars violate the condition that E{p) be monotonically decreasing. We expect this 
violation because the colormagnetic term for light pseudoscalars gives a large positive contribution 
to dE/dfi. 

Using the semi-empirical mass formula, we obtain a splitting in the B c system of 65 ± 10 MeV, 
and in the bb system of 55 ± 10 MeV. We estimate that the masses of the B c and r]b are 

M(B C ) = 6255 ± 30 MeV, M{r, b ) = 9405 ± 15 MeV. (35) 

We give the mass of the B c in Table II. 

We next turn to excited vector meson states, for which the data are considerably poorer 
than for the ground states. Furthermore, complications might arise from possible mixing with 
four-quark, hybrid, and glueball states. For example, the excited p(1465) and w(1394) states 
differ in mass by about 70 MeV, although they ought to be degenerate in our model. We believe 
this difference indicates appreciable mixing. Similar considerations apply to light P-wave mesons. 
Therefore, we use only the charmonium and bottomonium data of the Particle Data Group [22] 
and confine ourselves to excited B* states. Because we have only two data points for each excited 
state (a ip and an T) we cannot do better than use a linear fit to predict the masses of missing 
states, again using the quark masses in the first row of Table I. The problem with a linear fit is 
that a straight line is not concave upward, so that the predictions made in this fashion should be 
regarded as upper limits. We show in column 4 of Table II our predictions for the upper limits of 
two excited vector meson B* states. 

Quigg [36] suggested that it might be better to interpolate between mass differences, since 
these are considerably smaller than the masses themselves. With this procedure we have no 
theoretical reason to reject a linear interpolation. We show in column 5 of Table II our predictions 
for the masses of vector B* excited states using linear interpolation of the mass differences between 
corresponding states in the cc and bb systems. Again the errors include not only statistical errors 
but an estimate of the errors associated with the procedures. Note that the predicted masses in 
column 5 are less than the upper limits of column 4, i.e., the latter are indeed an upper bound. 

Turning to the P-wave mesons, we are not able to show analytically that the sum of the 
Fermi-Breit tensor and spin-orbit interactions satisfies the inequality (5). Nevertheless, it turns 
out empirically that if we use the same values of quark masses as for the vectors (row 1 of Table I), 
the eigenenergies of the tensor (J p = 2 + ), axial vector (J p = 1+), and scalar (J p = + ) mesons 
satisfy (9) and (25). We exploit this fact to fit separate 3-parameter exponential curves to the 
tensors, axial vectors, and scalars so as to obtain predictions for the B c P-wave states. These are 
shown in column 4 of Table II. Our estimated errors are rather large because of deviations of the 
curve from the eigenenergies of the known mesons. 

We can also use linear interpolation between cc and bb states to obtain the masses of B c 
P-wave states. These are also given in column 5 of Table II together with estimated errors. We 
see from Table II that the predictions for the P-wave mesons in columns 4 and 5 agree within the 
errors. 
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VIII. BARYON MASSES 

Because the masses of only four baryons of spin 3/2 are known experimentally and none 
of these contains any heavy quarks, the Feynman-Hcllmann theorem by itself does not enable 
us to make useful predictions of the masses of any baryons containing heavy quarks. However, 
if we use the Feynman-Hcllmann theorem in conjunction with a semi-empirical formula for the 
colormagnctic splitting in baryons [15], we are able to make some useful estimates of unknown 
masses. The reason is that the masses of the A c (quark content qqc), S c (qqc), S c (qsc), and Af, 
(qqb) spin-1/2 baryons are known from experiment [22,37,38], so that we can estimate the masses of 
the corresponding spin-3/2 baryons from a semi-empirical formula for the colormagnctic splitting 
in baryons [15]. We are then able to use a procedure analogous to that we used for mesons in order 
to obtain estimates of the masses of unknown baryons. 

The expression (32) for mesons can be generalized to baryons as follows [15]: We order the 
quarks so that if two quarks have the same flavor, they are chosen to be the first two; if all quark 
flavors are different, then the first two are the lightest. We denote by M* the mass of the ground 
state spin-3/2 baryon, with Ms the mass of the ground state spin-1/2 baryon in which the first 
two quarks have spin 1, and with Ma the mass of the spin-1/2 baryon whose first two quarks are 
in a relative spin state. We take the two-quark colormagnetic matrix elements as 

8Ri jt k = F i:j . k pas{2^ij)^ q i:j / {nii + mj), (36) 

with 

Fij,k = [Hij + x(mk + Mjfe)]/(M»j + A*ifc + Hjk) (37) 

to simulate the shrinking of the wave function with increasing mass nik of the spectator quark. 
The expression for in Eq. (37) is slightly different from that given in Ref. [15] and fits the 
data somewhat better. Following [15], we write the ground-state baryon mass differences as: 

M*-Afs=3fli3,2 + 3i223,i. 

M s - M A =4fli2, 3 - 2^13,2 - 2fl 23 ,i. 

The parameters n/ and Aqcd are chosen to be the same as in the meson case, while p, q, and x 
are adjustable. 

Although the structure of (38) is motivated by perturbation theory, we take these semi- 
empirical mass formulas to be more generally applicable, the justification being the good agreement 
with observed baryon mass splittings. 

A difficulty is that, unlike in the meson case, where we determined the quark masses from a 
fit to the vector meson eigenenergies, we do not know a priori what input values to use for quark 
masses which will be best suited for baryons. Our procedure is to start with m q = 300 McV and 
the other quark masses taken at reasonable starting values, for example, with the values given in 
[15]. We then adjust the quark masses by an iteration procedure which we shall now describe. 

(1) We use input quark masses and adjust the parameters p, q, and x to get a best fit to the 
known colormagnetic splittings in baryons. We then use the semi- empirical mass formula to 
calculate the masses of three spin-3/2 baryons (£*, S*, and ££) which are not known from 
experiment. 

(2) We then use the Feynman-Hcllmann theorem, analogously to the meson case; i.e., we adjust 
the parameters of a three-parameter curve and the quark masses m s , m c and nib to get a best 
fit to the eigenenergies of the A, £*, S*, 0, £*, £*, and ££. 

(3) We then use the new quark masses in the semi-empirical mass formula and repeat steps (1) 
and (2). 

In practice, this method rapidly converges. We find the best parameters of the semi-empirical 
mass formula (38) are p = 0.331 (GeV) 2 " 9 , q = 0.417, x = 3.805 when used with the following 
(rounded) quark masses for baryons in MeV: 



rriq = 300, m s = 475, m c = 1640, m b = 4990. 



(39) 
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The mass difference m s — m q satisfies the inequality (21). The parameters of an exponential curve 
of form (26) turn out to be 

a = 1307 MeV, b = 757 MeV, c = 813 MeV. (40) 

There are several reasons why our procedure for baryons is not as precise as that we used for 
mesons. First, the input "data" for baryons include three baryon masses which do not come from 
experiment but are only estimated from a mass formula. Second, even using these three baryons, 
we have only seven baryons to obtain the three parameters of a curve and the three quark masses. 
Third, we have to obtain unknown baryon masses by extrapolation, which is a less precise method 
than the interpolation method used for mesons. 

Comparing the masses of Eq. (39) with those used for mesons (see the first row of Table I), 
we see that the quark masses which give a best fit to the baryons are (except for m q , which was 
assumed to be the same) a little higher than those which lead to a best fit to the mesons. Because 
these quark masses are constituent masses, i.e. effective ones, there are no theoretical reasons why 
the masses determined from the baryons should coincide exactly with those determined from the 
mesons. If we insist that a single set of quark masses hold for both baryons and mesons, and vary 
these masses, our overall best fit to the hadron data is significantly poorer and our predictions 
have greater errors. 

Using the baryon data only, we can predict the masses of as yet unobserved baryons from the 
Fcynman Hcllmann theorem and the baryon semi-empirical mass formulas. As we have remarked, 
the baryon masses, given in Table III, are obtained by extrapolation, rather than interpolation, so 
that the errors are larger than in the meson case. The errors in the masses of the E b {qsb, spin 
1/2, antisymmetric in qs), E' b (spin 1/2, symmetric in qs), (spin 3/2), £l c (ssc, spin 1/2), fl* 
(spin 3/2), Cl b (ssb, spin 1/2), and f2£ (spin 3/2) arise partly because of the substantial error in 
the measurements to date [37,38] of the mass of the A;,. We believe that the following predicted 
mass differences (in MeV) are likely to have smaller errors than any of the masses given in Table 
III. 

M(£ 6 ) - M(A 6 ) = 200 ± 20, M(££) - M(A b ) = 230 ± 20, 

M(E b ) - M(Ab) = 190 ± 30, M(E' h ) - M(A h ) = 330 ± 30, (41) 

M(E* b ) - M(A b ) = 360 ±30. 

It is interesting that our semi-empirical mass formula makes the T, b about 10 MeV heavier than 
the Sfc. However, the probable error is such that this is not a firm prediction. All we can really 
say is that the T, b and E b have masses which are very likely within 20 MeV of each other. 

IX. CONCLUSIONS 

In conclusion, we have shown that the Fcynman Hcllmann theorem leads to the inequality 
SE/dfi < for most ground-state vector mesons and spin 3/2 baryons, even in the presence of 
relativistic kinematics and a flavor-dependent colormagnetic interaction. We were not able to 
show this for hadron pairs which differ by one member of the pair containing both a heavier and 
a lighter quark than the other, but the result seems to be empirically true even in this case. This 
inequality and the concavity condition (25), provide theoretical justification for an interpolation 
method [9], which allows one to make a quantitative prediction about the mass of the B* and 
other mesons without assuming any specific functional form for the quark-antiquark interaction. 
We obtain the masses of still other mesons containing heavy quarks by using a semi-empirical 
mass formula and by interpolating among mass differences. For the baryons, we can also use (8) 
and (25) to obtain predictions, but we need the semi-empirical mass formula from the outset and 
also need to extrapolate, rather than interpolate, in order to obtain useful results. Therefore, our 
predictions are not as precise as in the meson case. In making our predictions of the values of 
heavy hadron masses, we have not had to assume an explicit form for the Hamiltonian, but only 
some general characteristics about its flavor dependence. 
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TABLE I. Values of quark constituent masses in MeV for calculating meson energy eigenvalues 
from experimental values of their masses. We show in the first row the (rounded off) values of 
the quark masses used in this work and, for comparison, values used by some other authors in 
subsequent rows. 



Reference 


m q 


m s 


m c 


m b 


This work 


300 


440 


1590 


4920 


[9]* 


263 


404 


1543 


4876 


[15] 


300 


500 


1800 


5200 


[18] 


220 


419 


1628 


4977 


[21] 


310 


620 


1910 


5270 


[25] 


336 


510 


1680 


5000 


[26] 


337 


600 


1870 


5259 


[7] 


350 


500 


1500 


4700 


[14]* 


270 


600 


1700 


5000 


[17] 


335 


450 


1840 


5170 


[27] 


330 


550 


1650 


4715 


[28] 


150 


366 


1320 


4749 



*One of several sets of quark masses in this reference. 



TABLE II. Predicted masses of as yet unobserved B c (be) mesons. In column 4 we show 
predictions for the ground-state vector and P-wave mesons and upper limits for two excited vectors 
from interpolation of the energy eigenvalues, using the Feynman-Hellmann theorem. In column 5 
we show the pseudoscalar mass obtained from a semi-empirical mass formula, and excited vector 
and P-wave states from interpolation of mass differences. 



Name Spin-parity J p n 2S+1 Lj Mass (MeV) Mass (MeV) 



B* 


1" 


1 3 Sj 


6320 ±20 




B c 


o- 


1 % 




6255 ± 30 


B* 


1" 


2 3 Si 


< 6940 


6900 ± 20 


B* 


1" 


3 3 Si 


< 7290 


7250 ± 20 


B* 


0+ 


1 3 Po 


6630 ± 40 


6660 ± 30 


B* c 


1+ 


1 3 P! 


6730 ± 40 


6740 ± 30 


B* 


2+ 


1 3 P 2 


6770 ± 40 


6780 ± 30 
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TABLE III. Predicted masses of as yet unobserved baryons. In column 3 we show a prediction 
for a ground-state spin-1/2 baryon whose first two quarks have an antisymmetric spin wave 
function. (Antisymmetric O c and fib states do not exist in our picture.) Column 4 shows predictions 
for ground-state spin-1/2 baryons with symmetric spin wave function in the first two quarks. In 
column 5 we show predictions for the ground-state spin-3/2 baryons. See Eq. (41) for (we believe, 
more reliable) estimates for baryon mass differences. 



t Input 



Name 


Quark content 


M A (MeV) 


M s (MeV) 


M* (MeV) 


At, Eh, 


qqb 


t5630 ± 40 


5830 ± 40 


5860 ± 40 


77 77' 77* 


qsb 


5820 ± 40 


5960 ± 40 


5990 ± 40 


fl c ,fl* c 


ssc 




2710 ±50 


2770 ± 50 




ssb 




6070 ± 60 


6100 ± 60 



Figure captions 

FIG. 1. Energy eigenvalues of vector mesons using experimental masses from Refs. [22-24] in 
conjunction with four sets of quark masses from Table I: (a) from Ref. [21], (b) from Rcf. [18], (c) 
from Ref. [7], and (d) from Rcf. [27]. In the figure, the letters a through i stand for p, K* , <fi, D* , 
B* , £>*, B*, J/if), and T respectively. 

FIG. 2. Energy eigenvalues of vector mesons using our quark masses from the first row of 
Table I and the experimental masses from Refs. [22-24]. The letters stand for the same mesons 
as in Fig. 1, and the solid circle is our prediction for the B*. The solid line is a fit to the vector 
meson data with an exponential form, Eq. (26), with parameters given in Eq. (29). 

FIG. 3. Eigcncncrgies of the pseudoscalar mesons obtained from the masses of the Particle 
Data Group [22] (crosses) compared with the eigenenergies from the vector meson masses and the 
semi-empirical mass formula of Eq. 32) (open circles). The quark masses of the first row of Table 
I were used to obtain eigenergies from masses. In order of increasing /x are ir, K, D, B, D Sl B s , 
rj c , B c , and r\ h . 
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